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We determine the 9 dependence of the deconfinement temperature of SU(3) pure gauge theory, 
finding that it decreases in presence of a topological 6 term. We do that by performing lattice 
simulations at imaginary 8, then exploiting analytic continuation. We also give an estimate of such 
dependence in the limit of a large number of colors TV, and compare it with our numerical results. 
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The possible effects of a CP violating term in Quantum 
ChromoDynamics (QCD) have been studied since long. 
Such term enters the Euclidean lagrangian as follows: 



QCD 



iOq(x) 



■ 

?(*) = J^e^ pa F* v {x)F a pa {x) (1) 

where q(x) is the topological charge density. 

Experimental upper bounds on 9 are quite stringent 
(1^1 ~ 10~ 10 ), suggesting that such term may be forbid- 
den by some mechanism. Nevertheless, the dependence 
of QCD and of SU(N) gauge theories on 9 is of great 
theoretical and phenomenological interest. 9 derivatives 
of the vacuum free energy, computed at 9 — 0, enter vari- 
ous aspects of hadron phenomenology; an example is the 
topological susceptibility x = (Q 2 ) /V {Q = J d 4 x q(x) 
and V is the space-time volume) which enters the so- 
lution of the so-called U(1)a problem [U 0]. Moreover 
it has been proposed [H that topological charge fluctu- 
ations may play an important role at finite temperature 
T, especially around the deconfinement transition, where 
local effective variations of 9 may be detectable as event 
by event P and CP violations in heavy ion collisions. 

In the present work we study the effect of a non-zero 
9 on the critical deconfining temperature T c , considering 
the case of pure Yang-Mills theories. Due to the sym- 
metry under CP at 9 = 0, the critical temperature T c {9) 
is expected, similarly to the free energy, to be an even 
function of 9. Therefore we parameterize T c (9) as follows 



T c (0) 
To(0) 



1-R e 9 2 + Q{9 4 ) 



(2) 



In the following we shall determine Re for the SU(3) 
pure gauge theory, obtaining Rg > 0, and compare it 
with a simple model computation valid in the large TV 
limit, showing that Rg is expected to be 0(1/TV 2 ). 

The method - Effects related to the topological 9 term 
are typically of non-perturbative nature, hence numerical 



simulations on a lattice represent the ideal tool to explore 
them. However, it is well known that the Euclidean path 
integral representation of the partition function 



Z(T,6) = J [dA] e 



-S QC n[A]+i0Q[A\ _ p -V s f(8)/T 



(3) 



is not suitable for Monte- Carlo simulations because the 
measure is complex when 9^0. Sqcd — J d 4 x Cqcd 
and periodic boundary conditions are assumed over the 
compactified time dimension of extension 1/T; f(9) is 
the free energy density and V s is the spatial volume. 

A similar sign problem is met for QCD at finite baryon 
chemical potential /is, where the fermion determinant 
becomes complex. In that possible partial solu- 

tion is to study the theory at imaginary fiB, where the 
sign problem disappears, and then make use of analytic 
continuation to infer the dependence at real /is , at least 
for small values of hb/T [4]. An analogous approach has 
been proposed for exploring a non-zero 9 [5H8|]; as for 
fXB 7^ 0, also in this case one assumes that the theory is 
analytic around 9 = 0, a fact supported by our present 
knowledge about free energy derivatives at 9 = 0, ■ 

Various studies have shown that the dependence of the 
critical temperature on the baryon chemical potential, 
T c ((i>b)} can be determined reliably up to the quadratic 
order in \xb, while ambiguities related to the procedure of 
analytic continuation may affect higher order terms . 
It is natural to assume that a similar scenario takes place 
for analytic continuation from an imaginary 9 = i9i term, 
i.e. that Rg can be determined reliably from numerical 
studies of the lattice partition function: 



Z L (T,9) = J \dC] 



,-S L [U]-6 L Q L [U] 



(4) 



where [dU] is the integration over the elementary gauge 
link variables U^; Sl and Ql are the lattice discretiza- 
tions of respectively the pure gauge action and the topo- 
logical charge, Ql = J2 X We will consider the Wil- 
son action, Sl = PY^x,n>u(^ ~ ReTr II M ^ (x) /TV) where 
= 2TV/<7o an< 4 is the plaquette operator. 
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Various choices are possible for the lattice operator 
qL,{x), which in general are linked to the continuum q(x) 
by a finite multiplicative renormalization [l2j 



q L (x) a ^a 4 Z(P)q(x)+0(a 6 ) 



(5) 



where a = a((3) is the lattice spacing and lim a ^.o Z = 1. 
Hence, as the continuum limit is approached, the imagi- 
nary part of 9 is related to the lattice parameter 9l ap- 
pearing in Eq. (0} as follows: 9j = Z9^. 

Since q^ (x) enters directly the functional integral mea- 
sure, it is important, in order to keep the Monte-Carlo 
algorithm efficient enough, to choose a simple definition, 
even if the associated renormalization is large. Therefore, 
following Ref. j8|, we adopt the gluonic definition 



Ql(x) 



±4 



-1 



2 9 7T 2 ^ , 



Ti(U^(x)U pa (x)) , (6) 



where e 



^pvpa 



for positive directions and e 



pupa 



—£{-v)vp<j- With this choice gauge links still appear lin- 
early in the modified action, hence a standard heat-bath 
algorithm over SU (2) subgroups, combined with over- 
relaxation, can be implemented. 

Finite temperature SU (N) pure gauge theories possess 
the so-called center symmetry, corresponding to a multi- 
plication of all parallel transports at a fixed time by an el- 
ement of the center Zm- Such symmetry is spontaneously 
broken at the deconfinement transition and the Polyakov 
loop is a suitable order parameter. Since qb{x) is a sum 
over closed local loops, the modified action Sl + GlQl is 
also center symmetric, hence we still expect Zjq sponta- 
neous breaking and we will adopt the Polyakov loop and 
its susceptibility as probes for deconfinement 

Nt 



XL = 



V s ^ N 

x t=l 

V s ((L 2 ) (L) 2 )) 



(7) 



where Nt is the number of sites in the temporal direction. 

Results - In the following we present results obtained 
on three different lattices, 16 3 x 4, 24 3 x 6 and 32 3 x 8, 
corresponding, around T c , to equal spatial volumes (in 
physical units) and three different lattice spacings a ~ 
1/(4T C ), a ~ 1/(6T C ) and a ~ 1/(8T C ). That will permit 
us to extrapolate Rg to the continuum limit. 

We have performed, on each lattice, several series of 
simulations at fixed 6l and variable j3. Typical statistics 
have been of 10 5 — 10 6 measurements, each separated 
by a cycle of 4 over-relaxation + 1 heat-bath sweeps, for 
each run; autocorrelation lengths have gone up to O(10 3 ) 
cycles around the transition. In Fig. [1] we show results 
for the Polyakov loop modulus and its susceptibility as a 
function of /3 for a few values of 9l on a 24 3 x 6 lattice; 
we also show data obtained after reweighting in f3. We 
notice a slight increase in the height of the susceptibility 
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TABLE I: Collection of results obtained for j3 c and T c . 
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FIG. 1: Polyakov loop and its susceptibility as a function of ft 
on a 24 3 x 6 lattice and for a few 9l values. The susceptibility 
values have been multiplied by a factor 250. 



peak as 9l increases, however any conclusion regarding 
the influence of 9 on the strength of the transition would 
require a finite size scaling analysis and is left to future 
studies. 

The critical coupling /3 c (9l) is located at the max- 
imum of the susceptibility through a Lorentzian fit 
to unrewcighted data: values obtained at 9^ = 
coincides within errors with those found in previous 
works [13]. From /3 C (9 L ) we reconstruct T C (0 L )/T C (O) = 
a(/3 c (O))/a(/3 c (0i)) by means of the non-perturbative de- 
termination of a(/3) reported in Ref. [13(. Notice that 
most finite size effects in the determination of (3 c (9l) 
should cancel when computing the ratio T c (9l)/T c (0). 
A complete set of results is reported in Table HI 

As a final step, we need to convert 9l into the physical 



3 



parameter 9 = 19 j. A well known method for a non- 
perturbative determination of the renormalization con- 
stant Z = Z{fi) is that based on heating techniques fl4j . 
Here we follow the method proposed in Ref. [|[ , giving Z 
in terms of averages over the thermal ensemble: 



Z = 



(QQi 
(Q 2 ) 



(8) 



where Q is, configuration by configuration, the integer 
closest to the topological charge obtained after cooling. 
Such method assumes, as usual, that UV fluctuations 
responsible for renormalization are independent of the 
topological background. Z has been determined for a 
set of /3 values on a symmetric 16 4 lattice, as reported 
in Fig. [21 then obtaining Z at the critical values of f3 
by a cubic interpolation. Typical statistics have been of 
10 5 measurements, each separated by 5 cycles of 4 over- 
relaxation + 1 heat-bath sweeps, for each f3; the auto- 
correlation length of Q has reached a maximum of 10 3 
cycles at the highest value of (3. A check for system- 
atic effects has been done by repeating the determina- 
tion with a different number of cooling sweeps to obtain 
Q (15, 30, 45 and 60) or, at the highest explored value 
of (3, on a larger 24 4 lattice. In this way we finally ob- 
tain di(/3 c (6 L )) = Z(/3 C (9 L ))9 L , as reported in the 4th 
column of Table |U 

Final results for T c (9j)/T c (0) and for the three different 
lattices explored are reported in Fig. [31 In all cases a 
linear dependence in 9 2 , according to Eq. ([5]), nicely fits 
data. In particular we obtain Rg — 0.0299(7) for N t = 4 
(x 2 /d.o.f. ~ 0.3), Rg = 0.0235(5) for N t = 6 (xVd.o.f. ~ 
1.6) and Rg = 0.0204(5) for N t = 8 (x 2 /d.o.f. ~ 0.7). 

We have performed various tests to check the stability 
of our fits. If we change the fit range, e.g. by excluding, 
for each N t , the 1-2 largest values of 9i, results for Rg 
are stable within errors. If we assume a generic power 
like behavior T C (0)/T C (O) - 1 = A 9 a , we always obtain 
that a is compatible with 2 within errors; if we fix a to 
values which would imply a non-analyticity at 9 = 0, e.g. 
a = 1, we obtain a x 2 /d.o.f. of O(10) or larger. 

Assuming 0(a 2 ) corrections we can extrapolate the 
continuum value Rg = 0.0175(7), x 2 /d.o.f. ~ 0.97 (see 
Fig.[4j. Our result is therefore that T c decreases in pres- 
ence of a real non-zero 9 parameter. This is in agree- 
ment with the large N expectation that we discuss in the 
following, as well as with arguments based on the semi- 
classical approximation discussed in Ref. [IH for N = 2 
and with model computations [la ]. 

Large N estimate - We present now a simple argument 
to estimate the dependence of T c on 9 in the large N 
limit. Since the transition is first order, around the crit- 
ical temperature we can define two different free energy 
densities, f c (T) and fd(T), corresponding to the two dif- 
ferent phases, confined and deconfined, which cross each 
other at T c with two different slopes. The slope difference 
is related to the latent heat. Indeed the energy density 
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FIG. 2: Determinations of the renormalization constant Z on 
a 16 4 lattice. The dashed line is a cubic interpolation of data. 
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FIG. 3: T c (0)/r c (O) as a function of 6 2 for different values of 
Nt. Dashed lines are the result of linear fits, as reported in 
the text, then extrapolated to 8 2 > 0. 



is 



t 2 d 

( = V s df logZ] z = cxp 



Vsf(T) 
T 



(9) 



hence e = —T 2 d(f /T)/dT. Close enough to a first order 
transition we may assume, apart from constant terms, 
f c /T = Ad + 0(t 2 ) and f d /T = A d t + 0{t 2 ), where 
t = (T — T c )/T c is the reduced temperature. The latent 
heat is therefore Ae = — e c = T C (A C — Ad). 

A non-zero 9 modifies the free energy, at the lowest 
order, as follows: 

f(T, 9) = f(T, 9 = 0)+ x(T) 9 2 /2 + 0(9 4 ) (10) 

where x(T) is the topological susceptibility. x(T) is in 
general different in the two phases, dropping at deconfinc- 
ment [l7l - [l9j , hence the condition for free energy equilib- 
rium, f c = fd, which gives the value of T c , will change 
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FIG. 4: Re as a function of 1 /N? . The point at 1/Nt = is the 
continuum limit extrapolation, assuming 0(a 2 ) corrections. 



as a function of 9. The dependence of \ on T simplifies 
in the large N limit, being independent of T in the con- 
fined phase and vanishing in the deconfined one [HI, Eil . 
Hence we can write, for N — » oo, 



f c /T~A c t+( x /T) 9 2 /2 



fi/T ~ 4rf (11) 



where x is, from now on, the T = topological suscepti- 
bility. The equilibrium condition then reads (A c — Ad) t = 
( X /T c )9 2 /2 + 0(9 i ), giving 



Tc{6) 



= 1 - 



X 



O(0 4 ) 



(12) 



T c (0) * 2Ae 
In the large N limit we have 0, El Hflj], 

~ 0.0221(14) ; ~ 0.344(72) ; ~ 0.5970(38) 

apart from 1/A 2 corrections, hence we get 



Re 



X 
2Ae 



0.253(56) 
A 2 



0(1/N 4 



(13) 



The leading 1/N estimate for SU(3) is then Rg ~ 
0.0281(62). This is larger than our determination, even if 
marginally compatible with it: a possible interpretation 
is that for SU(3) the behavior of x a t T c is smoother 
than the sharp drop to zero that we have assumed. 

Notice that the 1/N 2 dependence of Re is in agreement 
with general arguments 21] predicting the free energy 
to be a function of the variable 6/N as N — > oo (see 
also Refs. [HEll)- For the same reason we expect 0{9 4 ) 
corrections to Eq. (fl2|) to be of 0{1/N A ): they are indeed 
related to 0(9 4 ) corrections to the free energy, which have 
been measured at T = by lattice simulations 22- 24j 
and are known to be small and of order 1 /N 2 . 

It would be interesting to extend the present study 
to N > 3, in order to check the prediction in Eq. (TT3"1) . 



and to N 
Ref. El- 



2, in order to compare with the results of 



We conclude with a few remarks and speculations re- 
garding the phase structure in the T~9 2 plane. In Fig. [3] 
we have drawn the critical line, for different N t and up 
to 6 2 terms, as fitted from 9 2 < simulations, and its 
continuation to 9 2 > 0; however other transition lines 
may be present, as it happens for the T — /i 2 ^ plane. 
For fi% < one finds unphysical transitions, known as 
Roberge- Weiss lines [25[, which are linked to the period- 
icity of the theory in terms of imaginary \i& ■ In the case 
of a 9 parameter, no periodicity is expected for imagi- 
nary 9, CP inviariance being explicitely broken for any 
9 1 7^ 0, hence we cannot predict other possible transi- 
tions for 9 2 < 0. A 27r periodicity is instead expected 
for real values of 9, with the possible presence of a phase 
transition at 9 = it where CP breaks spontaneously. 

Our simulations have given evidence, for 9 2 < 0, only 
for a deconfinement transition line, describable by a 9 2 
behavior up to \9\ ~ tt. We expect continuity of such be- 
havior, at least for small real 9, while non-trivial correc- 
tions may appear as 9 approaches 7T. However, following 



Ref. [2l[ and the arguments above, we speculate that, 



at least for large N, T c (9) be a multibranched function, 
dominated by the quadratic term down to 9 = n 



T C (0)/T C (O) 



1 - R e min (9 + 2nky 



(14) 



where k is a relative integer: in this case periodicity in 
9 implies cusps for T c {9) at 9 = (2k + 1)tt, where the 
deconfinement line could meet the CP breaking transition 
present also at T — 0. Therefore the phase diagram 
at real 9 could have some analogies with that found at 
imaginary fiB- 
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